Black holes in a box: 
Towards the numerical evolution of black holes in AdS space-times 
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The evolution of black holes in "confining boxes" is interesting for a number of reasons, particularly 
because it mimics the global structure of Anti-de Sitter geometries. These are nonglobally hyperbolic 
space-times and the Cauchy problem may only be well defined if the initial data is supplemented 
by boundary conditions at the timelike conformal boundary. Here, we explore the active role that 
boundary conditions play in the evolution of a bulk black hole system, by imprisoning a black hole 
binary in a box with mirrorlike boundary conditions. We are able to follow the post-merger dynamics 
for up to two reflections off the boundary of the gravitational radiation produced in the merger. 
We estimate that about 15% of the radiation energy is absorbed by the black hole per interaction, 
whereas transfer of angular momentum from the radiation to the black hole is only observed in 
the first interaction. We discuss the possible role of superradiant scattering for this result. Unlike 
the studies with outgoing boundary conditions, both the Newman-Penrose scalars *I'4 and ^^o are 
nontrivial in our setup, and we show that the numerical data verifies the expected relations between 
them. 

PACS numbers: 04.25.D-, 04.25.dc, 04.25.dg, 04.50.-h, 04.50.Gh, 04.60.Cf, 04.70.-s 



I. INTRODUCTION 

Numerical relativity and the gauge/gravity correspon- 
dence are two fields in high energy/gravitational physics 
that have seen tremendous activity and progress over 
the last few years. Since the 2005 breakthroughs [iHll 
the numerical relativity community has produced stable 
evolutions of black hole (BH) binaries in 4-dimensional, 
asymptotically flat space-times, covering the inspiral, 
merger and ring-down phases. State of the art simu- 
lations can be as long as 15 orbits plus merger and ring- 
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down [4], can deal with large mass ratios {^-Q, with di- 
mensionless BH spins of up to 0.92 Q or with eccentric- 
ities as low as 5 X 10^^ (see [l3| for a recent review). 
Such simulations have unveiled new phenomena such as 
BH kicks [ll|-[l§| and some groups are now focusing on 
the construction of sufficiently accurate template wave- 
forms to be used in filtering data from the current (LIGO, 
Virgo, GEO, TAMA) and planned (Advanced LIGO and 
LISA) gravitational wave observatories Simul- 
taneously. high energy collisions of BHs have been simu- 
lated [18l - l20| . leading to accurate results for the scatter- 
ing cross section and for the conversion rate of the initial 
centre-of-mass energy into gravitational radiation. These 
simulations have also tested cosmic censorship, the Dyson 
luminosity limit and exhibited the zoom-whirl behaviour 
first found in [21}. 

The gauge/gravity correspondence has been devel- 
oped, since the original proposal by Maldacena in 1998 
[23 | , both as an efficient tool to understand strongly cou- 
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pled gauge theories using classical gravity and as way 
to study gravitational phenomena from a dual field the- 
ory. Central to many of these developments are BHs, as 
may be seen by the following list of examples: (i) the 
successful microscopic computations of the Bekenstein- 
Hawking entropy for extremal BHs 12311 and Hawking 
emission rates for near extremal BHs |24| are now seen 
as applications of the correspondence; (ii) the confine- 
ment/deconfinement phase transition in QCD-like theo- 
ries has been identified [l^ with the Hawking-Page phase 
transition for Anti-de-Sitter (AdS) BHs jla]; (iii) mov- 
ing away from thermal equilibrium, the quasinormal fre- 
quencies of AdS BHs have been identified with the poles 
of retarded correlators describing the relaxation back to 
equilibrium of a perturbed dual field theory [l^, ^\ ; (iv) 
in a large class of gauge theories with a gravity dual, 
a universal behaviour was obtained for the ratio of the 
strongly coupled medium's viscosity to entropy density, 
by computing the absorption cross-section of low en ergy 
gravitons in the dual BH (or black brane) geometry [23] ■ 
The result is in good agreement with experimental re- 
sults from the Relativistic Heavy Ion Collider (RHIC); 
(v) critical exponents, of the type found in spherical grav- 
itational collapse by Choptuik in 4-dimensional asymp- 
totically flat space-time [30|, have been conjectured to 
be dual to the asymptotic value of the parton satura- 
tion exponent for high energy scattering in QCD, in the 
Regge limit, at weak coupling [3l|. These and other ex- 
amples have built expectations that the gauge/gravity 
correspondence will be a useful laboratory for gaining un- 
derstanding of very difficult problems of both field theory 
and gravity, such as the confinement problem, the infor- 
mation loss paradox and the problem of singularities. 

Given the potential of the correspondence exemplified 
above, a working framework to solve Einstein's equations 
exactly for a broad range of initial conditions in AdS 
spaces would, unquestionably, be very useful for deep- 
ening the study of the correspondence, particularly in 
more dynamical situations. To assemble such a frame- 
work is our long-term goal and the present paper serves 
as the first step in this direction, i.e., to generalize the 
techniques of numerical relativity to AdS space-times. 
In order to achieve this goal, one has to go beyond the 
standard methods of numerical relativity in, at least, two 
obvious points. 

Firstly, AdS space-times are not globally hyperbolic. 
In asymptotically AdS spaces the boundary plays an "ac- 
tive role" for the bulk evolution. This is easily visualised 
in the Penrose diagram of AdS, which has a timelike 
boundary. Physically, null geodesies in AdS reach the 
boundary for a finite affine parameter. One thus often 
refers to an asymptotically AdS space as a "box", hav- 
ing in mind that AdS boundary conditions directly affect 
the bulk physics f32| - [34| . This should be contrasted with 
the asymptotically flat case, where the only physically 
relevant choice for the boundary conditions of the bulk 
fluctuations corresponds to outgoing waves at spatial in- 
flnity. In the gauge/gravity correspondence, the choice of 



the AdS boundary conditions is dictated by a holographic 
prescription |35l - l38| . Secondly, from the viewpoint of the 
duality, I?-dimensional AdS space-times, AdS/j (not just 
AdS4) are relevant. Thus, we would like to have a frame- 
work that could be used in AdS d , h^ particular for D — 5, 
which is related by the correspondence to 4-dimensional 
gauge theories. The latter issue has been recently ad- 
dressed by our group in separate publications 39, 40] , as 
well as by other groups |4l| - |46| . Here we shall focus on 
the former issue: the active role of boundary conditions. 

The dynamics of BHs in AdS, and especially the role of 
spatial infinity ( "the box" ) is poorly understood. In con- 
trast to the asymptotically flat case, interesting new phe- 
nomena may occur in AdS backgrounds. For instance, 
superradiance effects have been shown to make small 
(as measured by the AdS radius) rotating BHs unstable, 
through a sequence of reflections at the boundary and 
amplifications close to the ergoregion [47l - l55| . The final 
state of this instability could be a new nonaxisymmet- 
ric BH configuration, which is also suppor ted by recent 
gravity/hydrodynamics arguments [56j, |57[. Notice that 
nonaxisymmetric BHs are strictly forbidden in asymp- 
totically flat space-times [H, [11] , so the boundary does 
have an important role in the description of BHs. 

In order to identify in the cleanest possible way the 
active role of the boundary for the bulk evolution, we 
consider here a toy model for AdS. We set the cosmo- 
logical constant to zero and impose mirrorlike boundary 
conditions on a box that contains the dynamical system. 
This mimics the AdS global geometry, keeping the local 
geometry of vacuum models. We choose the dynamical 
system to be a BH binary, starting at some given dis- 
tance, producing either a head-on collision or an inspi- 
ralling merger. In the latter situation we consider the 
initial BHs without intrinsic angular momentum. These 
are, by now, very well tested systems when purely out- 
going boundary conditions are imposed. Thus we will be 
able to see clearly the modiflcations due to the nonout- 
going boundary conditions in systems with nontrivial dy- 
namics. 

For the post-merger dynamics the inspiralling binaries 
provide a more interesting analysis than the head-on col- 
lision case, since the initial centre-of-mass energy trans- 
ferred into gravitational radiation is by more than one 
order of magnitude larger in the former case as compared 
to the latter. Immediately after the merger, the system 
will contain a single (spinning or nonspinning) BH plus 
gravitational radiation. This radiation will then be (re- 
peatedly) reflected off the boundary and interact with 
the BH. 

The flrst nontrivial result is that we can follow the 
numerical evolution for up to two reflections off the 
boundary of the gravitational radiation produced in the 
merger. A priori it was not guaranteed that this could 
be achieved, since it is not known whether the formula- 
tion of the Einstein equations that we use provides a well 
defined initial value boundary problem together with the 
boundary conditions we impose. Our simple setting ac- 
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tually provides a first attempt to test the well-posedness 
of the initial boundary value problem in a non globally 
hyperbolic space-time. We find that our numerical re- 
sults are at least second-order convergent for at least two 
reflections off the wall, after which we gradually lose con- 
vergence. A deeper study of these issues is clearly needed, 
as well as an exploration of how the convergence (and 
remaining results) change for different boundary condi- 
tions. 

During the window of numerical convergence, we study 
the properties of both outgoing and ingoing gravitational 
radiation. The usual studies of BH binaries with outgo- 
ing boundary conditions, focus only on the Weyl scalar 
^1/4, which describes outgoing gravitational waves. How- 
ever, an equally relevant quantity for the description of 
gravitational radiation is the scalar ^I'o, which describes 
ingoing waves, but which is seldom discussed in the liter- 
ature. Due to our special boundary conditions and setup, 
we are able to verify certain relations between these two 
quantities in a numerical evolution for the first time. This 
also provides a test on the correctness and meaning of the 
boundary conditions we have imposed. 

By analysing the properties of the apparent horizon of 
the BH produced in the merger and after each interac- 
tion with the gravitational wave packet, we estimate the 
amount of energy and angular momentum that is trans- 
ferred from the radiation into the BH per interaction. In 
case of the inspiralling binary the boxed BH is spinning 
and we expect superradiant scattering of the waves gen- 
erated during merger to become important; in fact, the 
back and forth bouncing of the waves at the reflecting 
wall and their subsequent amplification by superradiance 
close to the ergoregion are expected to turn the system 
into a BH bomb [H, [H, leoljH Thus, these simulations 
will be the first attempt at a nonlinear study of the BH 
bomb. An important open problem is understanding how 
the evolution proceeds and what is the end point of the 
instability. This can only be achieved through nonlinear 
studies. A final statement on this issue will, however, 
require further analysis than that provided herein. 

This paper is organized as follows. Section [IT] briefly 
reviews the setup to evolve Einstein's equations numer- 
ically for the case at hand, including a brief description 
of the numerical code, of the formulation of Einstein's 
equations in a so-called BSSN form, of the gauge choice 
adopted, the way the "spherical" boundary is imposed 
and the boundary conditions. In Section Hill we explain 
which of the numerical outputs we use to extract rele- 
vant physical quantities, in particular gravitational wave 



The artificial mirror sometimes appears naturally. A massive 
scalar field scattering off a Kerr BH acts as its own refiecting 
wall [6ll467tl . Furthermore, Kaluza- Klein modes in dimensional 
reduction can also act as an effective mass rendering higher- 
dimensional, rotating BHs unstable I68l l69ll . Finally, it has been 
suggested that astrophysical BHs might sometimes behave as BH 
bombs, with the role of the refiecting cavity being played by ac- 
cretion disks [70l[7l|. 



estimates (Section IIII A[) and apparent horizon (AH) es- 
timates (Section IIII Bp . The numerical results of our 
simulations are shown in Section IIVI In Section |V] we 
close with some discussion of the results and prospects 
for the future. Some technical points and further re- 
sults have been organised into three appendices. For self- 
containedness. Appendix |A] reviews the electromagnetic 
decomposition of the Weyl tensor and in particular the 
construction of the relevant quantities for our study, 
and ^"4. Appendix iBl exhibits some snapshots for visu- 
alising the evolution of the system we have studied. Ap- 
pendix [C] describes the simulations with a cubic, rather 
than spherical, box. 



II. NUMERICAL FRAMEWORK 

In order to numerically generate a solution to the Ein- 
stein field equations, it is most convenient to view the 
problem as a time evolution or initial value problem. The 
majority of formulations of the Einstein equations as an 
evolution system in time is based on the canonical "3-1-1"- 
decomposition introduced by Arnowitt, Deser and Mis- 
ner (ADM) ^ and further developed by York flB\. One 
thus obtains a first order, constrained evolution system 
in time for six components each of the three-metric jij 
and the extrinsic curvature Kij which describe the intrin- 
sic geometry of three-dimensional hypersurfaces as well 
as their embedding in the four-dimensional space-time. 
The Hamiltonian and momentum constraints impose four 
conditions on 7^ and Kij on each hypersurface but are 
conserved under the time evolution. Finally, four gauge 
variables, the lapse a and the shift represent the co- 
ordinate freedom of Einstein's relativity. Suitable speci- 
fication of these free variables is crucial for a successful 
numerical implementation. 

Our numerical framework is based on a method now 
commonly referred to as moving punctures 0, 0]. The 
Einstein equations are formulated as the Baurngarte- 
Shapiro-Shibata-Nakamura (BSSN) system [tI, uM' ^ 
modification of the ADM formulation which employs the 
variables 

X = '4'^* = 7"* , 7»j = X7y , 
K = Y^K.j , = xA^J = X (k,j - p^JK^ , 

r = 7^"'=f;.fe = -d,r ■ (1) 

By construction det jij — 1 which implies the last equal- 
ity for . The exact form of the evolution equations for 
this set of variables is given in Eqs. (Al, A4, A6, A7, 
Asflof Ref. [Z^. Finally, we evolve the gauge variables 



^ Note that the final term on the right hand side of their Eq. (A6) 
should be |x(<^^ ~ 9m/3™), i- e. a factor of x is missing. 
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a and /?' using "1+log" slicing and a F-driver of the form 



dta = (3"'drna - 2aK, 



(2) 
(3) 
(4) 



Here Xl3 s-^d VP ^tre constant parameters set to one for 
all simulations reported in this work. 

We evolve these equations with the Lean code 76 1 
which is based on the Cactus computational toolkit 77 1 
and the Carpet mesh refinement package [zl, [t^] . BH 
binary initial data are provided by the spectral solver of 
Ansorg et al. [s^ and the calculation of apparent hori- 
zons is performed with Thornburg's AHFinderDirect 
[Sll |82| . For more details of the code we refer the reader 
to Ref. [zil and Sec. Ill of [83^. 

The key ingredient in which our current numerical 
framework differs from previous implementations of the 
Lean code and most other codes is the outer boundary 
condition, which we will discuss in more detail in the 
remainder of this section. 

The vast majority of numerical simulations of BH bi- 
naries has been concerned with asymptotically flat space- 
times and consequently employed either of the following 
boundary treatments: (i) outgoing Sommerfeld condi- 
tions on Cartesian grids of finite size, as described for ex- 
ample in [i^l ; (ii) outgoing radiation [H, [11] with multi- 
patch methods, including Cauchy characteristic wave ex- 
traction [13, HI] and (iii) constraint preserving boun dary 
conditions combined with multidomain methods (89ll90j|. 

In contrast we will study the dynamics of BH space- 
times under the influence of a reflective outer boundary. 
It is natural to use for this purpose an outer boundary 
of spherical shape. Most importantly, this avoids mix- 
ing of different gravitational wave multipoles as would 
occur in the case of a reflective, cubic outer boundary. 
This is discussed in more detail in Appendix [C] where 
we compare simulations using both types of boundary. 
Except for this comparison, however, we will exclusively 
study spherical outer boundaries or, rather, approximate 
these by using so-called Lego spheres; cf. Sec. 3 in [oH . 
In Fig. [2] we sketch the foliation of the space-time un- 
der consideration, suppressing one spatial dimension for 
simplicity. The numerical domaine, i.e., the Lego sphere 
is visualized by a dark (red) domain on each timeslice 
St+„5t. Their numerical implementation is illustrated in 
Fig. [1] which schematically displays a computational do- 
main using four refinement levels with one or two compo- 
nents each. The individual components are labelled GJ^ 
where the indices i and m denote the refinement level 
and component number. Note that one spatial dimen- 
sion is suppressed for visualisation purposes. In order 
to update a grid function at a particular vertex, we re- 
quire information from neighbouring points because of 
the discretization of spatial derivatives in the evolution 
equations. The exact number n of neighbouring points 
required in each direction depends on the finite difference 
stencils employed. While n = 3 for the 4*^^ order accurate 



stencils used in our simulations, we use n = 1 for simplic- 
ity in our illustration in Fig. [TJ Consider first the dark 
(blue) shaded area inside the inner solid circle of radius 
Rb- Each point in this regular domain can be updated 
straightforwardly provided we also have valid data on 
the boundary points marked by x symbols. Points out- 
side the circle of radius Rb are not required for updating 
regular points and are simply ignored in the numerical 
evolution. The specific boundary condition is then deter- 
mined by the manner in which we update grid functions 
on the boundary points marked as x in the figure. 

In order to mimic the global structure of an Anti-de 
Sitter space-time we effectively enclose the BH binary 
inside a spherical mirror and set 



d_ 

dt 



/ = 0, 



(5) 



at each boundary point with / denoting any of the BSSN 
variables listed in Eq. ([IJ . The use of fourth-order sten- 
cils adds one complication to this picture: the upgrade of 
a grid point requires two neighbors, so that points right 
next to the boundary need special treatment. In practice, 
we have achieved optimal stability properties by evolving 
these points with second-order stencils. Our implemen- 
tation requires one further ingredient in order to handle 
the spurious radiation inherent to numerically generated 
initial data of BH binary systems; cf. [l^]- In order to 
avoid contamination of our simulations by such spurious 
radiation being trapped inside our refiective boundary we 
employ standard outgoing radiation boundary conditions 
at early times and only switch on our reflective condition 
at 



tret = i?B + At 



pulse ■ 



(6) 



In order to avoid a discontinuous jump from outgoing to 
reflective boundary conditions, we gradually switch off 
the time derivative df /dt using a weighting factor w{t) 
which smoothly decreases from 1 to over an interval 
At = 10 M and At = 20 M for the head-on collision and 
inspiral, respectively. The duration of the spurious wave 
pulse Atpuise is estimated from previous simulations of 
similar setups in asymptotically flat space-times as for 
example presented in Refs. [tgI, [ol, [OJ]. The spurious 
radiation is thus given sufficient time to leave the com- 
putational domain. 



III. EXTRACTION OF MEANINGFUL 
QUANTITIES 

A. Wave extraction 

Information about gravitational radiation is most com- 
monly obtained in numerical simulations of BH binaries 
by computing the Weyl scalar ^'4, which has the asymp- 
totic property of being equal to the outgoing radiation 
if the complex null-tetrad is chosen properly. In a truly 
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FIG. 1. Illustration of a (Lego-)spherical outer boundary. 




FIG. 2. Sketch of the foliation for the numerical evolution of 
BH binaries in a (spherical) box. The location of the consid- 
ered numerical domain on each spatial hypersurface is shown 
as a dark (red) sphere. 



To be explicit, we define a spherical coordinate sys- 
tem centred on the centre of mass of the binary with 
orthonormal basis {f,d,(j)). The coordinates are chosen 
such that the azimuthal axis is aligned with the orbital 
angular momentum and the binary orbits are in the direc- 
tion of increasing azimuthal coordinate. Our definitions 
and notation are the same as in [tB, [9^ . To define our 
complex null-tetrad, we use the timelike unit vector nor- 
mal to a given hypersurface h and the radial unit vector 
f to define an ingoing (fc) and outgoing null vector {£) 
by 

k = -^(h + f), i=-^(h-f). (7) 

We define the complex null vector m and its complex 
conjugate by 

m=^{$ + i9), m=^{4>^ie). (8) 

In terms of this tetrad, we define and ^'4 as 

4-0 = Co.0^sk"m^k''m\ (9) 
*4 = Ca/s^^rm'^FTO^ (10) 

where Cap^s is the Weyl tensor. To relate ^'o and 
to the amplitudes of the gravitational waves, we note 
that in the transverse-traceless (TT) gauge, assuming the 
functional form f{t±r) for the ingoing or outgoing waves, 
we have 



1 



R^ 



TT 



-hi 



R^ 



±R^ 



2 ■ "9(1) "ridnip " rOrtp nOf(f) 



±R^ 



A 'dot' denotes derivative with respect to the argument. 
Following standard conventions, we take the /i+ and hx 
polarisations of the gravitational waves to be given by 



TT 



TT 



(11) 



asymptotically AdS space-time this procedure cannot be 
implemented so straightforwardly. In our toy model, 
which has the same local geometry as vacuum models, for 
a sufficiently large box size we expect the standard tetrad 
to yield the correct gravitational wave information. This 
is yet another advantage of our setup. Moreover, because 
we also deal with ingoing waves in our simulations, once 
they are reflected from the box boundary, we will con- 
sider the Weyl scalar V^o as well, in order to account for 
the ingoing contribution. The authors are not aware of 
any numerical study of the ingoing Weyl scalar ^fp. We 
will therefore analyse it in detail in the course of our 
numerical studies. 



Then, we find that in vacuum regions of the space-time, 
for outgoing waves /i+ = /i+(t — r) and ~ hx{t ^ r), 



*o = 0, 
4-4 = /m 



(12) 
(13) 



while for ingoing waves = h^[t+r) and = hx (t+r) 



^0 = /^^ 

^4 = 0. 



ihy 



(14) 
(15) 



The fact that 4*4 (4*0) are zero for ingoing (outgoing) 
waves is consistent with the leading order in perturba- 
tion theory. The correct expression and in particular its 
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dependence on the radial coordinate is given by the solu- 
tion of the Teukolsky master equation at large distances 
from the source, which states that for outgoing waves [96l | 



while for ingoing waves 

*0 « 



„iw(t+r) 



(16) 



(17) 



We decompose the resulting ^'4 {'^0) into modes by pro- 
jection onto spherical harmonics of spin- weight s = — 2 
(s = 2) according to 

00 / 

MrvI/4 = Mr ^ ^ _2i^™(^,^)V'L, (18) 

OQ I 



1=2 m=-l 



where -2Yim{0, (p) and 2Yim{9, cj)) are spin- weight —2 and 
2 spherical harmonics [97, 98]. These are defined as 



(20) 



where d\^g is the Wigner c?-function 



^ {-iy^{l + m)\{l-m)\{l + s)\{l-s)\ 



t=Ci 



{l + m-ty.{l 



t)W.{t 



my. 



(cos 61/2) 



2l+m-s-2t 



(sin 61/2) 



2t-\-s — 7n 



(21) 



and where Ci = max(0, m — s) and C2 = min(Z-|-m, l — s). 
Here M is the ADM mass of the system, computed from 
the initial data and assuming this is an asymptotically 
fiat space-time, and r is the generalised harmonic radial 
coordinate. 

In the numerical code, the null-tetrad is constructed 
from a Cartesian orthonormal triad (u, w, w) and the 
timelike vector h is orthonormal to t = constant hy- 
persurfaces. The space-time is evolved with time t using 
Cartesian coordinates x, y, z. In practice, we compute 
the Newman-Penrose scalars \l/o and ^4 using the elec- 
tromagnetic decomposition of the Weyl tensor accord- 
ing to Eqs. (|All|) on the entire Cartesian grid. Then, 
they are interpolated onto coordinate spheres of differ- 
ent extraction radii r^^ with a uniform distribution of 
points in (0, (/<) All the waveform related data from 
the simulations presented in the course of this paper are 



^ 108 X 54 points in S e [0, tt] , G [0, 27r] for the set of simula- 
tions using low resolution. In case of the higher resolutions these 
numbers are adjusted accordingly. 



taken from such samplings of ^'o(i, r = fex,^, 0) and 
^'4(t,r = rex, 0,0). A more detailed description is given 
in Appendix [Al 

As discussed previously, 5*4 is no longer simply related 
to the energy flux, but since we are dealing with fairly 
large box sizes, one might hope that many notions can be 
retained in an approximate sense. Given the Newman- 
Penrose scalar ^'4, we can compute the radiated energy, 
linear and angular momentum from the radiation con- 
tent |99|]: 



dE 
~dt 



—— = lim 



dt 

dJ^ 
dt 



16tt 



lim 

iim X 



^idi 



dn, 

2 



-^Adt 



dn. 



(22) 
(23) 



Re 



^Adt a 




-ilAdtdt dn 



— 00 ^ — 00 



(24) 



where 



£i — {— sin 6 cos cj), — sin sin (p, — cos ( 



The definitions above are based on time integrals which 
start in the infinite past (at retarded time t = —00), 
and thus capture the complete gravitational wave signal. 
Starting the time integrations at t = —00 corresponds 
to the limit of infinite extraction radius on the initial 
time slice — the slice would then extend all the way to 
spatial infinity, no part of the waveform would be lost, 
and it would take an infinite time for the waves to reach 
the extraction sphere. This situation cannot be handled 
with the current numerical codes; we therefore work with 
finite extraction radii. 

The mass and angular momentum of the final BH can 
be estimated from balance arguments. Given the param- 
eters Py- , d in the Bowen-York initial data, we straight- 
forwardly calculate the total initial angular momentum 
as 



Jir 



dP„. 



(25) 



since the initial spin of each BH is zero. Ansorg's TwoP- 
UNCTURES [iOl initial data solver directly provides the 
total ADM mass M of the system and we obtain radi- 
ated energy and angular momentum -Erad and Jj-ad from 
the gravitational wave signal. In case of a merger, this 
gives us the final angular momentum and mass of the BH 



Jtin 



(26) 
(27) 



The dimensionless spin parameter of the final hole follows 
directly from 



Jfin 



Jfin 



(28) 
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We check our results by fitting the quasinormal frequency 
and damping time of the final BH and invert them to 

obtain j'qnm (see e.g. [innl - [in2} v 



B. Apparent horizon properties 

We can also characterise the process by the proper- 
ties of the apparent horizon of the final BH itself. Since 
this relies only on local quantities, it does not depend 
upon the space-time being asymptotically flat. In order 
to monitor the mass and spin of the final BH, we use 
Thornbur g's Apparent Horizon Finder AHFinderDi- 
RECT [sil |82| in different ways, which also allow us to 
obtain uncertainty estimates: 

1. The irreducible mass Mi„ enables us to calculate 
the final BH mass Mbh from Christodoulou's rela- 
tion [Toi 



J' 



(29) 



This relation provides a method to check the inter- 
nal consistency of the result for the final BH spin as 
calculated from the above balance arguments. For 
this purpose we set A/bh = -^^fin and solve Eq. (p9)) 
for the spin 



Jfin 



ML 



= 4 



ML 



ML 



(30) 



For comparison we also compute the spin of the 
final hole from the two following estimates: 

2. We measure the ratio Cr{j) = Gp/Ce of polar to 
equatorial circumference of the final BH [104] . If we 
assume the final object to be a Kerr BH, this ratio 
is Cr = l^l- P'^EiP^), where = fM/{2r+), 
E{0^) is a complete elliptic integral and r+/M — 
1 



p. This expression can be inverted to 



find the dimensionless spin parameter, j^^, of the 
final hole. 

The equatorial circumference of a Kerr BH is C, 
47rM. Therefore 2t: Aah/CI = I + ^1 
where Aah is the area of the apparent horizon. 
Thus, the AH area and the equatorial circumfer- 
ence can be u sed to estimate the spin of the final 
BH from [Tol 



Jah^ 



JAH = -1/1- 



2t:A 



AH 



- 1 



(31) 



IV. NUMERICAL RESULTS 

Our numerical study focuses on two types of binary 
BH initial configurations; (i) head-on collisions of non- 
spinning BHs starting from rest and (ii) quasicircular in- 
spiral of nonspinning holes. In the remainder of this work 
we label these as HD and IN simulations. The initial pa- 
rameters of all our simulations as well as the structure of 
the computational domain and the position of the outer 
boundary Rb are summarised in Table HI Unless denoted 
otherwise, the results presented refer to the highest res- 
olution available. 



A. Numerical convergence analysis 

Before we discuss in detail the physical properties of 
the BH binary systems, we test the accuracy of our sim- 
ulations by performing a convergence analysis of model 
IN2. Numerical simulations based on the finite dif- 
ferencing method typically approximate the continuum 
solution of differential equations with a leading error 
that has a polynomial dependence on the resolution, 
/cont = /num + O (/»"). The Order n depends on the spe- 
cific numerical implementation. Consistency of the code 
is tested by evolving the same configurations with low, 
medium and high resolution /ic, /im and h / . One straight- 
forwardly shows that the convergence factor is then given 
by 



_ /he — !h 



fh„, - fhf 



(32) 



where fh stands for any of the evolved variables ob- 
tained for resolution h. For the case of contemporary 
moving puncture codes, the solution is complicated by 
the fact that the differential equations are typically dis- 
cretized with fourth- (or higher) order accurate stencils 
but prolongation in time between different refinement 
levels and implementation of outer boundary conditions 
is only second-order accurate; see Sec. IV in Ref. |l07l | 
and our discussion in Sec. [Hi 

Fig. [3] shows our convergence results for the vari- 
ables ^'4 and \l/o obtained for resolutions he = M/48, 
h,n = M/52 and hf = M/56. Here the differences 
fh,„ — fhf have been amplified by the convergence factors 
Q — 1.26 and 1.47 expected for second and fourth-order 
convergence, respectively. The figure indicates fourth- 
order accuracy for the first passage of the wave pulse and 
a gradual deterioration of convergence to second-order 
accuracy afterwards. We believe this to be a consequence 
of the different ingredients of the code as discussed above. 
At early stages, the dominant error is the discretization 
of derivatives. As the pulse successively passes across 
mesh refinement boundaries and is reflected off the outer 
boundary, however, the second-order error in the pro- 
longation operation becomes dominant and reduces the 
order of convergence. We also note, in this context, that 
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Run 


Grid Setup 


Rb/M 


d/M 


Mirr,i/M 


Pr/M 






INI 


{(48, 24, 12, 6) X (1.5, 0.75), h = 1/56} 


40 


6.517 


0.483 


±0.133 


0.69 


0.70 


IN2.1 
IN2.2 
IN2.3 


{(48, 24, 12, 6) X (1.5, 0.75), h = 1/48} 
{(48, 24, 12, 6) X (1.5, 0.75), h = 1/52} 
{(48, 24, 12, 6) X (1.5, 0.75), h = 1/56} 


30 
30 
30 


6.517 
6.517 
6.517 


0.483 
0.483 
0.483 


±0.133 
±0.133 
±0.133 


0.69 


0.65 
0.65 
0.65 


HDl 


{(48, 24, 12, 6) X (1.5, 0.75), h = 1/60} 


40 


6.517 


0.483 


0.0 


0.0 


0.0 


VIS 


{(48, 24, 12, 6) X (1.5, 0.75), h = 1/48} 


48 


6.517 


0.483 


±0.133 







TABLE I. Grid structure, as well as initial and final parameters of the simulated black holes. The grid setup is given in terms 
of the radii of the individual refinement levels as well as the resolution near the punctures h (see Sec. II E in [7y] for details). 
The reflective outer boundary is located at radius Rb- The tab le further shows the initial coordinate separation of the two 
punctures d, the irreducible mass Mirr and the Bowen-York [l06l | parameter for initial linear momentum Pi of the individual 
holes. and Jgn"^ are the spin of the single hole after merger determined from the AH and the merger radiation, respectively. 

All parameters are given in units of the ADM mass. We did not monitor the AH properties for all runs, therefore the spin of 
the final BH is not determined (empty cells) for some cases. The Weyl scalars have been extracted at rex = 35M (INI, HDl) 
and Tex = 25M (IN2), respectively. Model VIS is used in Appendix [B] for visualisation. 




FIG. 3. Convergence analysis of the outgoing Weyl scalar ^4 (left panel) and the ingoing Weyl scalar ^0 (right panel) for 
the IN2 runs. We show the differences of the I = m = 2 mode between the coarse and medium and the medium and fine 
resolution run. The latter has been amplified by the factors Q = 1.47 (fourth-order convergence) and Q = 1.26 (second-order 
convergence). We observe fourth-order convergence in the signal due to the merger whereas the first and second after- merger 
cycles show only second-order convergence. The first two refiected and ingoing wave pulses show second-order convergence. 



well-posedness of the BSSN evolution system with re- 
flective boundary condition has so far not been demon- 
stratecfl. We can therefore not rule out adverse effects 
on the long-term convergence properties due to potential 
ill-posedness of the continuum system of equations. In 
the remainder of this discussion we will restrict ourselves 
to 2-3 passages of the wave pulse as covered in Fig. [3] 



To our knowledge, the well-posedness of the system of equations 
in combination with reflecting boundary conditions, as treated 
here, has not been studied yet. Some investigations of the wave 
equation with th is type of boundary conditions suggest that it 
may be ill-posed llOSl Il09ll . These investigations also show that 
the wave equation with periodic bou ndary co nditions is a well- 
posed initial boundary value problem Il08l . ll09l . pointing towards 
interesting future extension of our work. 



during which the relative uncertainties in ^4 and ^0 are 
< 5%. 



B. Gravitational wave signal and black hole 
dynamics 

To our knowledge, this work presents the first analy- 
sis of gravitational waveforms with both outgoing (^'4) 
and ingoing (^'o) contributions for long-term stable nu- 
merical simulations of BH binaries. For this reason, we 
first illustrate the general pattern of the wave signal ob- 
tained for model VIS of Table HI A series of snapshots 
of both Newman-Penrose scalars are shown in Fig. IH] in 
Appendix |B] in superposed form. 

The gravitational wave signal is dominated by the 
quadrupole contributions which we show in Fig. 21 For 
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100 200 300 400 

t/M 



FIG. 4. Real part of the / = m = 2 mode of rM>I>o and rM*4 
of run INI. The ingoing signal tM'^q has been shifted in time 
by At = lOM and in phase by tt (thus equivalent to an extra 
minus sign) to account for the additional propagation time 
and the reflection. 



clarity, the ingoing signal ?A22 has been shifted in time 
by At = 10 Af in order to compensate for the additional 
propagation time from the extraction radius rex = 35 M 
to the boundary Rb = 40 M and back after reflection. 
The reflection introduces an additional phase shift of 
A0 — TT which has also been taken into account in the 
figure. Within numerical errors, we find the resulting 
outgoing and subsequent ingoing pulses to overlap. 

The first outgoing wave pulse, visible in Fig. |4] around 
150 < t/M < 250, is generated during the inspiral, 
plunge and merger of the binary and is similar to wave- 
forms obtained for the inspiral of nonspinning BH bi- 
naries in asym ptotically fiat space-times (cf. Fig. 1 in 
Refs. [tI, Il07l |). Due to the reflecting boundary, how- 
ever, this wave pulse does not escape the computational 
domain. Instead it propagates inwards, interacts with the 
post-merger remnant hole and eventually manifests itself 
as a second wave pulse shifted by At w 80 M relative 
to the first. This process repeats itself many times, with 
the wave pulse being presumably distorted (by absorp- 
tion, superradiance and other curved space-time effects 
on wave propagation) upon each interaction with the BH. 
We now investigate in detail these changes of the wave 
pulse upon interaction with the BH. 



1. Interaction of the wave pulse with the remnant black hole 

As shown in Fig. HI the outgoing and subsequent ingo- 
ing wave pulses overlap within numerical uncertainties. 
We therefore focus on the outgoing signal in our study of 
subsequent wave pulses and the gradual changes caused 
by successive scattering off the BH. Changes in the wave 
pulse are best illustrated by considering the wave ampli- 
tude as shown in Fig. [S) Here we superpose the 1 = 2, 



m = mode for model HDl and the I = 2, m = 2 
multipoles for models INI and IN2 of the first three suc- 
cessive outgoing wave pulses by applying corresponding 
time shifts to the waveform. Clearly, the wave pulses 
broaden after each scattering off the BH. We emphasise 
that this distortion of the pulse is not an artifact of the 
outer boundary condition as is demonstrated by the good 
overlap between the ingoing and outgoing pulses in Fig.U 

One possible explanation for this distortion relies on 
strong field, curved space-time effects: massless waves in 
curved geometries do not propagate on the light-cone. 
Indeed, as shown in the classical work by DeWitt and 
Brehme, the Green's function for a massless field in a 
curved spa ce-tim e does not generally vanish inside the 
light cone [llOl |. This gives rise to interesting effects. 
It is well-known that wave propagation in BH space- 
times shows that the signal can roughly be divided in 
three parts: (i) the first part is the prompt response, at 
very early times, whose form depends strongly on the 
initial conditions. This is the most intuitive phase, be- 
ing a counterpart of the light cone propagation in fiat 
space-time; (ii) at intermediate times the signal is dom- 
inated by an exponentially decaying ringing phase, and 
corresponds to the e xcita tion of the BH's characteristic 
modes of vibration [102^; (iii) a late-tim e tail, usually 
a power law falloff of the field Ill2j |. Therefore, a 
variety of possible distortions are possible. 

Another possible explanation for the increasing width 
of the pulse arises in the context of the superradiance 
mechanism. It is well-known that the scattering of a 
wave pulse off a rotating Kerr BH will result in super- 
radiant scattering — amplification of the scattered wave 
packet via extraction of rotational energy from the BH — 
if the wave pulse satisfies 

uj < mfl , (33) 

where = jfin/(2r+) is the BH angular velocity [6^, 
Ill3l Ill4| . We note, however, that Fig. [5] also indicates 
a broadening of scattered wave pulses in the head-on 
case where no superradiance is expected. While the ob- 
served broadening would be compatible with superradi- 
ance, other effects appear to also influence the shape 
of the pulse and our observations do not conclusively 
demonstrate its presence. 

In order to investigate this dispersion in more detail, 
we plot in Fig. [6] the energy spectrum for the domi- 
nant I = 2, m = 2 mode. The vertical line in this 
figure denotes the threshold frequency mft 0.4/Af cor- 
responding to a final spin jfin = 0.69 as obtained for 
the post-merger hole for configuration IN2; cf. Sec. lIIIBl 
The figure demonstrates that the I = 2, m = 2 mode 
does contain contributions which would be subject to 
superradiance-induced amplification. These results then 
suggest that the low-frequency component of the pulse is 
amplified due to superradiance, while the high-frequency 
component is absorbed. A linear analysis of superradi- 
ance in the Kerr geometry [ll4| shows that superradi- 
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FIG. 5. Overlap of the amplitudes of successive pulses of the same waveform; / — 2,m — for the HDl run (left), I = m = 2 
for the INI (centre) and IN2.3 (right panel) simulations, obtained by time-shifting such that the maxima overlap. 




Mco 



FIG. 6. The energy spectrum for the I — m — 2 multipole of 
the outgoing scalar ^'4, for model IN2.3. This spectrum corre- 
sponds solely to the initial wave packet, i.e., the waveform has 
been truncated immediately before the first reflection off the 
boundary. The vertical line marks the threshold frequency 
for superradiance. 



ant effects are always small, unless the hole is rotating 
close to the extremal value. Thus, further studies, includ- 
ing larger spins of the post-merger hole, are necessary to 
comprehensively demonstrate superradiant wave ampli- 
fication. 



2. Black hole dynamics 

In contrast to the case of asymptotically flat space- 
times, a BH binary under the influence of a reflective 
shell does not settle down into a stationary configuration 
soon after merger. This is due to the repeated interac- 
tion with the wave pulse passing back and forth across 
the spatially finite space-time. The prolonged dynamical 
state of the system manifests itself prominently in the 
area Aah of the apparent horizon. In Fig. [7] we show the 
fractional deviation (Aah ~ ^o)/^o of th^ apparent hori- 
zon area from its value Aq immediately after merger, i.e. 
the first instance a common apparent horizon is found. 



0.03 



0.02 




FIG. 7. Time evolution of the area of the apparent horizon 
for the head-on and inspiral simulation HDl (solid curve) and 
INI (dashed curve). The area of the rotating BH increases 
at regular intervals corresponding to the propagation time of 
the pulse between the hole and the reflective boundary at 
Tex = 35 M. Due to the small amount of radiation generated 
during the plunge in the head-on case (HDl), the variation in 
the AH area is buried in numerical noise. 



As expected, the horizon area remains nearly constant for 
the duration of the first passage of the pulse to the reflec- 
tive shell and back. At « 80 M and 60 M, respectively, 
for models INI and IN2.10. The subsequent increase in 
Aah demonstrates that some fraction of the gravitational 
wave pulse energy is absorbed by the hole. It follows an- 
other period of approximately constant horizon area, a 
further increase upon the second scattering of the pulse 
and so on. Unfortunately, the radiation efficiency is al- 
most 2 orders of magnitude lower for head-on collisions 
(see Table II in Ref. [iHl), so that the increase in hori- 
zon area is buried in the numerical uncertainties. The 



^ Small oscillations in the horizon area are due to numerical un- 
certainties. 
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FIG. 8. Time evolution of tlie (relative) mass of the BH (solid 
line) computed by M = d/'iTi, the irreducible mass (dashed 
line) and the total spin J = jM^ (dashed-dotted line). 



head-on case serves as a useful comparison, however, as 
it demonstrates that the changes observed for the inspi- 
ralling configurations are significant relative to numerical 
uncertainties. 

The BH mass, defi ned in terms of the equatorial radius 
of the horizon Ce by jl05j | 



Ce 
477 



(34) 



shows a similar behaviour as the horizon area. In Fig. [8] 
we plot the fractional deviation {M—Mq)/Mo of the mass 
from its value immediately after merger together with the 
irreducible mass and the BH spin J of the hole obtained 
for model INI. The mass remains approximately con- 
stant until the pulse returns after its first reflection, then 
increases, remains constant during the second passage of 
the pulse and so on. In contrast, the spin shows a signifi- 
cant increase only during the first scattering of the pulse 
off the BH. 

We conclude that in each interaction with the gravita- 
tional radiation, the final BH mass increases. It is inter- 
esting to compare the increase in the horizon mass with 
the amount of gravitational wave energy radiated dur- 
ing the last stages of the inspiral, plunge and merger of 
a corresponding binary system in an asymptotically fiat 
space-time which is about 3.5 % of the total energy of the 
system [t^ [q^I • For the INI run, we estimate that about 
15% of the energy emitted during the merger is absorbed 
by the central spinning BH per interaction. Our results 
are consistent with total energy conservation. More- 
over, they are not incompatible with superradiant am- 
plification: typically, absorption of high-frequency waves 
is more effective than superradiant amplification of low- 
frequency waves, such tha t the net effect leads typically 
to absorption by the BH |114| . The prolonged nonsta- 
tionary character of the post-merger state is also demon- 
strated by the time dependence of the BH's final spin. 



Immediately after merger, we obtain jfin = 0.69 from 
Eq. pip in excellent agreement with corresponding sim- 
ulations using outgoing radiation boundary conditions 
[q^ I- Successive interaction with the reflected wave pulse, 
however, results in a small but significant increase in the 
BH's spin as shown in Fig.[S] As before, the first increase 
occurs about At = 80 M after merger, when the pulse 
has returned to the BH. We estimate the fractional in- 
crease in spin resulting from the first scattering at about 
5 %. For comparison, the total angular momentum ra- 
diated in the case of an asymptotically fiat space-time is 
reported as J'^^^/M^ = 0.246 (28 % of the initial orbital 
angular momentum of the system) in Table I of Ref. . 
It thus appears that a significant amount of angular mo- 
mentum remains in the form of gravitational waves. Due 
to numerical uncertainties it is not entirely clear whether 
later periods of interaction between pulse and hole result 
in a further transfer of angular momentum from the wave 
pulse to the hole or vice-versa. Our results indicate, how- 
ever, that the amount of angular momentum exchanged 
in subsequent interactions is significantly below 5 %. 



V. DISCUSSION 

The dynamics of BHs in generic space-times is a 
fascinating, yet extremely challenging problem. The 
gauge/gravity duality, however, strongly motivates us to 
solve dynamical problems with BHs on asymptotically 
AdS backgrounds. In this paper, we have studied a toy 
model that captures one of the fundamental features of 
such backgrounds: the active role played by the bound- 
ary conditions for the bulk evolution. 

We have mimicked the global structure of an AdS back- 
ground by introducing a refiecting wall at some radius. 
Within this cavity we evolved an inspiralling BH binary 
and a BH binary starting from rest at a certain initial 
distance. Of course, these are very specific initial con- 
figurations, and serve merely as tests for future, possibly 
more complex, situations. 

Perhaps the most important conclusion of the present 
work is that these simulations can be done and represent 
the first step to a full numerical evolution of BHs in AdS 
space-times. Indeed, as observed in Section IIV Al it is 
not known whether the BSSN evolution scheme together 
with reflecting boundary conditions is a well-posed initial 
value problem. Thus, the convergence we have exhibited, 
which holds up to two reflections off the boundary of the 
gravitational radiation produced in the merger, is the 
first of our results. Among the other results presented 
here we stress the following: 

1. For the first time, we were able to numerically study 
the scalar V^o, describing ingoing waves, and check 
in the numerical data the simple relations between 
*o and *4 (cf. Eqs. dH]), ^ and Fig.®. 

2. Our results are consistent with the intuitive expec- 
tations for a wave packet of radiation (generated 
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during inspiral plus merger) travelling back and 
forth between the mirrorlike wall and the BH: part 
of this radiation is absorbed when interacting with 
the BH (especially high-frequencies). We estimate 
that about 15% of the wave packet's energy is ab- 
sorbed by the BH per interaction, at least during 
the first cycles. 



3. The wave packet is clearly distorted upon interac- 
tion with the hole, which means some frequencies 
are absorbed more efficiently than others. In prin- 
ciple, we should observe a small but nonzero su- 
perradiance. Unfortunately, we have not obtained 
incontrovertible proof of nonlinear superradiance 
from our numerical data. Presumably, the system 
will become unstable after a sufficiently long time, 
since the radiation should be exponentially ampli- 
fied, once the high-frequency components had time 
to be completely absorbed by the BH. In future 
work we plan to investigate these instability stud- 
ies further by considering a highly spinning, final 
BH produced by the inspiral of spinning BHs. 

One issue that we have not explored in this paper is the 
potential influence of the box on the premerger dynam- 
ics. For sufficiently small size boxes, it is plausible that 
the radiation produced in the inspiral may be refiected 
off the boundary and interact with the binary before the 
merger. This might produce observable signatures in the 
premerger dynamics and even in the properties of the 
remnant black hole. Whereas we seem to observe some 
hints of this effect in our numerical data, a more exhaus- 
tive analysis is required to produce some precise state- 
ments. 

A future interesting extension of the present work is 
to repeat our analysis for periodic boundary conditions. 
Indeed, the wave equation in combination with peri- 
odic boundary conditions is a well-posed initial boundary 
value problem [108, 109]. Such system is thus more likely 
to yield longer stable and convergent numerical evolu- 
tions, allowing one to address a number of interesting 
effects. Some relevant questions are: for how long is the 
numerical evolution stable? Do the numerical instabil- 
ities set in before the system has had time to reach a 
stable configuration? What are the typical times needed 
to achieve this? 

Also in the context of numerical relativity in space- 
times with a cosmological constant, we plan on inves- 
tigating dynamical black hole space-times in de Sitter 
backgrounds, where many of the pro blems present in AdS 
do not exist (see also Refs. |ll5l4ll7 for some work along 
these lines). 

Finally, it is necessary to extend this work in the ob- 
vious fashion: implement the evolution of BHs in real 
AdS backgrounds. We hope the present work will help 
on achieving that. 
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Appendix A: Electromagnetic decomposition of the 
Weyl tensor 



Since in this work we analyse the Newman-Penrose 
scalar ^fp, besides \E'4, which is uncommon in numerical 
works, we collect in this appendix some useful results. 
Following the sign convention in (tgI [95| the Newman- 
Penrose scalars and ^4 are computed by eqs. © and 
([T0|. The vectors k,{.,m,ifi form a null-tetrad. Their 
inner products vanish except for 



kt=\=mm. 



(Al) 



In practice, the vectors of the null-tetrad are constructed 
from a Cartesian orthonormal basis (u, w) in the spa- 
tial hypersurface and the timelike orthonormal vector h 
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according to 



1 



V2 
1 



(A2) 



The orthonormal triad vectors are constructed via the 
Gram- Schmidt orthonormaUzation starting with 



v' = [xz,yz,-x^ - y'^], 



(A3) 
(A4) 
(A5) 



where Cjf. is the 3-dimensional Levi-Civita tensor. Next, 
we decompose the Weyl tensor in terms of its electric and 
magnetic parts [ost 



(A6) 



where = 7^^ + n^n^ and Capj = e/^i/Ap"'" l.'^J-pJ-P. 
The electric and magnetic part of the Weyl tensor are 
given by 



Bc.p = *Cp.Ap K nP. 



(A7) 
(A8) 



_L(^ denotes the projection operator onto the hypersur- 
face and * denotes the Hodge dual. By using the Gauss- 
Codazzi equations we express the electromagnetic com- 
ponents in terms of the "3-1-1" variables [zl] 



E- — R 



-/""'{K.jKki-K.kKji), (A9) 

(AlO) 



If we insert Eq. (jA6[) . the definition of the null-tetrad 
(|A2p and the expressions (|A9I) . (jAlOp into the definition 
of the Newman-Penrose scalars (l^. (|10p we obtain 



^-4 = -[Eki{v''v-w''w') - Bkiiv^w' +v'w'')] 



(All) 



In the numerical code we use these relations in order to 
calculate ^E'o and ^"4 on the entire Cartesian grid. Then, 
they are interpolated onto coordinate spheres of various 
extraction radii rex- The Newman-Penrose scalars ^"0 
and ^4 are decomposed into spin- weighted spherical har- 
monics sYirn according to 



In the first equation we have used the relation [97 



.Y, 



Im 



(-) 



m-\-s 



-.Yu 



Thus, in practice we implement 
terms of the spherical harmonics 
-2. 



-2Ylm 



(A12) 



and 'i/'fm only in 
with spin- weight 



Appendix B: Snapshots 



In Fig. [9] we illustrate the emission of the gravitational 
wave signal during the inspiral and merger and its evo- 
lution in the closed (confined) system containing a cen- 
tral, spinning BH. We display snapshots of the wave- 
forms by superposing (the real part of) 4'o and ^'4 as 
obtained for model VIS of Table HI We show a slice of 
the orbital plane with x,y = —ASM, ...,48M during an 
interval t/M = 150,..., 540. The difference in time be- 
tween the individual pictures is At/AI = 10. The series 
of snapshots starts in the late inspiral phase shortly be- 
fore the plunge and we see a strong gravitational wave 
signal that is emitted throughout the merger (first row 
and first two columns of the second row). This signal 
reaches the spherical boundary and is reflected back as 
can be seen in the final three panels of the second row 
and first two panels of the third row. Starting with the 
third snapshot in the third row we see a second pulse 
going outwards again after it has been scattered off the 
BH. This process is repeated several times and the series 
of snapshots ends with the fourth outgoing wave pulse. 
An animati on c onstructed from the numerical data can 
be found in IllE 




FIG. 9. Snapshots of 5R(^'4) superposed by 3?(»I'o). The snapshots show the evolution from t — 150M until t — 540M and 
a time interval of At = lOM. We show a slice of the orbital plane with both coordinates going from — 48M, ASM. 
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Appendix C: Cubic-shaped box 

We have also performed simulations using condition 
([5]) on a cubic outer boundary. These simulations repre- 
sent the inspiral of nonspinning BH binaries with initial 
separation d/M = 6.514 and initial linear momentum 
PyjM ~ ±0.133. The grid setup for these runs is listed 
in Table Hi] together with a reference model 01 which 
describes the inspiral of the same binary in an asymptot- 
ically flat space-time using outgoing radiation boundary 
conditions. Gravitational waves have been extracted at 



Run 


Grid Setup 




RbIM 


b.c. 


Ol 


(192, 96, 48, 24, 12, 6)(1.5, 0.75) 


1/40 


192 





Cl.l 


(24,12,6)(1.5, 0.75) 


1/40 


24 


R 


C1.2 


(24,12,6)(1.5, 0.75) 


1/44 


24 


R 


C1.3 


(24,12,6)(1.5, 0.75) 


1/48 


24 


R 



TABLE II. Parameters for a set of models evolved using a cu- 
bical boundary with reflective boundary condition ("R b.c") 
for models Cl.l, CI. 2 and CI. 3 and an outgoing ("O b.c") 
Sommerfeld condition for model 01. 



Tex = 20 M in the form of the Newman- Penrose scalar 

In Fig. [ini we compare the I ~ m — 2 mode of ^"4 
obtained by the evolution of models Cl.l and 01. 



0.08 
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FIG. 10. Comparison of the I — 2,m = 2 modes of ^4 ob- 
tained for models Cl.l (solid) and Ol (dashed curve). The 
expected range in time for subsequent wave pulses resulting 
from first and second reflections are indicated in the figure. 



In Fig. [TT] we present the convergence plot of the 
I = m = 2 mode of ^'4 obtained from evolutions of models 
Cl.l, CI. 2 and CI. 3. The difference between the medium 
and fine resolution result has been amplified by the fac- 
tor Q — 1.58 corresponding to fourth-order convergence. 
While the overall convergence is about fourth-order as in 
the case of a spherical shell, the cubical outer boundary 
introduces a substantial amount of numerical noise which 
is also demonstrated in Fig. [T^ which displays the 1 — 2, 

m — 2 and the I = A, m — A multipoles of ^4 for models 
0.03 1 1 1 1 1 1 1 1 1 1 1 1 — 



0.02 



r'refl. 




300 



50 



100 



1 50 200 
t/M 



250 



300 



350 



FIG. 11. Convergence analysis of the Z = m = 2 mode of ^4 
obtained for model CI of Table HT] 



CI. 3 and INI. We believe that this is at least partly a 
consequence of mode mixing in the case of the cubical 
boundary shell which is not well-suited for the geometric 
shape of the gravitational wave pulse. For this reason, we 
have exclusively used a spherical shell in the main part 
of this work. 
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